Shape invariant potentials in SUSY quantum mechanics and periodic orbit theory 
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We examine shape invariant potentials (excluding those that are obtained by scaling) in super- 
symmetric quantum mechanics from the stand-point of periodic orbit theory. An exact trace formula 
for the quantum spectra of such potentials is derived. Based on this result, and Barclay's functional 
relationship for such potentials, we present a new derivation of the result that the lowest order 
SWKB quantisation rule is exact. 
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In non-relativistic quantum mechanics certain poten- 
tials are amenable to exact analytic solution. For a subset 
of these soluble potentials, the energy spectrum may be 
expressed explicitly as an algebraic function of a single 
quantum number. Such potentials occur either in one 
space dimension, or are central potentials in higher di- 
mensions. For the latter, an effective potential in the ra- 
dial variable can be defined for each partial wave. Some 
examples of such potentials are Coulomb, harmonic os- 
cillator, Morse, Rosen-Morse, etc. 0. These potentials 
also have the property that the lowest order WKB quan- 
tisation rule, together with the appropriate Maslov in- 
dex (that may change from potential to potential 0) , 
leads to exact results. For central potentials, the Langer 
prescription 3] for the centrifugal barrier, together with 
half-integer quantisation, can also be employed Q. 

In supersymmetric (SUSY) quantum mechanics, these 
exactly solvable potentials are found to be translation- 
ally shape invariant ||. Combining SUSY and WKB, 
Comtet et al. @ found that the lowest order SWKB cal- 
culation needs neither the Maslov index nor the Langer 
correction to yield the exact result. The purpose of the 
present paper is to understand this result from the point 
of view of the periodic orbit theory (POT) rather 
than the higher order WKB corrections Q. Regarding 
the latter, we should point out a largely overlooked pa- 
per by Barclay 0, in which he showed that the higher 
order WKB terms converge in these potentials to yield 
an energy-independent correction. The latter may be ab- 
sorbed in the Maslov index of the lowest order term. For 
SWKB, all the higher order terms vanish. Although we 
do not make the WKB expansion, we arrive at the same 
result in a novel application of POT. 

We first set the notation by reviewing the relevant 
equations of SUSY QM. Consider a potential V(x;ai) 
of a single variable x, and a set of parameters denoted 



by a\. One defines a "super potential" 

h <f/ (x) 



W{x;ax) 



/ 2m 0o (x) 



where 4>q(x) is the ground state solution of the 
Schrodinger equation at energy Eq for the potential 
V(x,ai), and a prime denotes the spatial derivative. Let 
us define 

V 1 (x;a 1 ) = (V(x;a 1 )-E Q ) , (1) 
so that the ground state energy of the Hamiltonian 



H 



2m dx 2 



Vi{x;ai) 



lies at zero energy, i.e., E^ = 0. Then it is easy to show 
that 

Vi{x; oi) = W 2 (x; oi) - -^W'(x; oi) . 



The SUSY partner Hamiltonian H 2 has the potential 
V 2 (x;ai) , and has an energy spectrum identical to that 
of Hi, except for the absence of the zero-energy state. 
The ground state of H 2 , denoted by Eq ' coincides with 
the first excited state E± of Hi, and so on. The partner 
potential V>2(x;ai) is 



-.W'(x; ai 



V 2 {x;a 1 ) = W 2 {x;a 1 ) + 



Shape invariance in the partner potentials is defined by 
the relation 



V 2 {x;a 1 ) = Vi(x;a 2 ) + R(ai) , 



(2) 



where the new parameters a 2 are some function of ai, 
and the remainder R(ai) is independent of the variable 
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x. We restrict our consideration of shape invariance to 
those cases where a% and a\ are related by translation, 
a 2 = a i + a - It is then straight forward to show, by con- 
structing a hierarchy of Hamiltonians that the complete 
eigenvalue spectrum of Hi is given by Q 



Ei 1 



4 13 = o 



(3) 
(4) 



The RHS of the above may be expressed as a monotonic 
function fi(n) of the quantum number n, so that 



E^=h(n); /i(0) = 



(5) 



For the shape invariant potentials we consider here, /i(n) 
is an algebraic function. Using this property, we proceed 
to obtain an exact expression for the quantum density of 
states of Hi in the spirit of periodic orbit theory. This en- 
tails a division of the density of states into a smooth and 
an oscillating part as a function of a continuous classical 
variable E. To this end, we may write 

5(E - £W) = S(E - /i(n)) = S(n - Fi(E))F{(E), (6) 

where the relation E = fi (n) has been inverted to define 

n = Fi(E). (7) 

(When the mapping between n and E n is (1, 1) there is 
no difficulty in the inversion. But to define the deriva- 
tive F[{E) requires an extension from the discrete to the 
continuum, and that is not unique. In the case of shape 
invariant potentials we have an algebraic relation which 
solves that problem.) For the spectrum under consider- 
ation, /i(0) = implies the condition 



Fi(0) = 



(8) 



The quantum density of states gi(E) for the discrete 
spectrum of Hi is defined as 



gi(E) = J2d(n)S(E-EP) 



(9) 



n=0 



where d(n) is the degeneracy of states at E = E n . Writ- 
ing d(n) = d(Fi(E)) = D(E), and using Eq. ©, we 
obtain 

oo 

gi(E) = D{E)F[{E) £ S(n - ft (25)) . (10) 

n=0 

(D(E) = 1 for one-dimensional potentials). We now use 
the identity 

oo oo 

J25(n-x)= ^ . *>°* (11) 



to obtain the desired expression 0, ^| 



gi(E) = D(E)F[(E) 



2^cos[27rfcFi(£;)] 



fc=i 



(12) 



For a given ft (25), this is an exact expression for the 
quantum density of states Q±(E). It is in the form of a 
trace formula in POT || when Fi(E) (to within a 
dimensionlcss additive constant jf) is related to the ac- 
tion Si (E) of the primitive classical periodic orbit of the 
potential Vi(x) : 



Si(E) 



= f 1 (e) + t I , 



rX2 

Si(E) = 2V2^i / y^E~Vidx 



(13) 
(14) 



In the above, xi and xi are the classical turning points 
at which E = V\ (x) (for simplicity in notation, we write 
Vi(x,ai) — Vi(x)). The (h— independent constant) r\ 
may be determined by using Eq. (|13fl . and applying the 
condition given by Eq. © for E = 0. We then obt am 



5i(0) 



(15) 



We may prove Eq. (|13fl by noting that the (smooth) 
Thomas-Fermi density of states, given by the first term 
on the RHS of Eq. l|12f) . is the Laplace inverse of the clas- 
sical canonical partition function |l2| of the Hamiltonian 

Hf(x,p) =p 2 /2m + Vi(x) : 

F{(E) = C- E l Zfifi) = — / Z?{0)eP E dp. (16) 



Since 



Zf((3) = lj exp[-pH?(x,p)] dxdp 



1 I2rrnr 



2nh V 13 
it follows from (|16|l that 



exp[-/3Vi(x)]dx , (17) 



F{(E) 



'2m 



dx 



2nh J X1 y/[E-Vi(x)] 



(18) 



From this, Eq. i|13|) follows on integration over energy. 
Using Eq. J7J) together with l|13l I14|) , we obtain the im- 
portant result that the lowest order WKB quantisation 
rule is exact for Vi : 



Si(E) = <b p{x)dx = (n + rj)h 



(19) 



n=0 



k= — oo 



where p(x) = \j2m\E — Vi(x)\. We also see that the 
constant r\ is the so called Maslov index which may vary 
from one potential to the other. 
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The Maslov index r\ may be eliminated from the quan- 
tisation rule by employing the superpotential formalism, 
and the result of Barclay and Maxwell . They made 
the important observation that the shape invariant class 
of potentials under consideration obey one or other of the 
following equations: 
Class I 

'' " W = A + BW 2 {x) + CW(x) , (20) 



\f2m dx 
or Class II 
h dW 
V2m dx 

where A, B and C are constants. Using these equations, 
we now show that Si(E), as defined by Eq. (|14|l . obeys 
the relation (xi s , X2 S are the turning points in SWKB) 

S^E) = 2^2^i / \JE-W 2 dx + hr) 



= A+BW 2 (x) + CW(x) y/(A + BW 2 ) , (21) 



(22) 



To this end, note that the action S\ can be expressed as 
an inverse Laplace transform 

-f3[W 2 -hW /y/2m\ 

Si(E) 



v2w7r £ E 1 



£3/2 



-dx . (23) 



At this point, for simplicity of notation, let us temporar- 
ily put %/ 1 \j2vn = 7. Expanding the exponential in pow- 
ers of W, we have 



v2m7r£. 



£ 

k=0 



nk+i 



(Sfe + iji 



2V2m / \JE - W 2 dx 



£ 



d k 



VE 



l\k 



{k + l)\8E k VE - w 2 



dW 



(25) 



Note that now the limits in x are replaced by the con- 
dition W 2 (x) = E. The integral for k = may be done 
immediately, yielding tt. To evaluate the integrals for in- 
teger k > 1, we assume that -fW' obeys Barclay's equa- 
tion <H3) (class I) or (J3TJ (class II). 

For class I, we require integrals of the type 



h.= 



(a + bw 2 + cwy 



dW 



(26) 



-Ve VE-W 2 
On expanding the numerator, terms with odd powers of 
W vanish on integration. One now sees that only the 
piece of Ik involving the highest power of W 2 survives 
the differentiation in Eq. I|25|) . Consider the integral 
with W 2k . With the substitution W = VI sin 



Ve 



W 



2 k 



-Ve VE - W 2 



-.dW = E k 



it/2 



sin 2k 6>d<9 



-tt/2 



^•(2fc-l)!! 



(2k)\ 



(27) 
(28) 



Accordingly, Eq. (|2*5|l reduces to 
Si(E) 



2v 2m / Ve - w 2 dx 

kl 



fc=l 



B k {2k- 1)!! 
(fe + l)(2Jfe)!! 



(29) 



By construction, W 2 (x) has coincident turning points at 
E = 0, so the first term on the RHS above vanishes at 
this energy. Comparing with Eq. JTSJl, we deduce that 



k=l 



B k (2k-1)\\ 
(fc + l)(2fc)!! 



(30) 



Note, from Eq. (|20|l . that B is independent of Planck's 
constant h. Comparing now with Eq. (|13|l . we deduce 
our main result 



2jrhF 1 {E) = V2m I \J E - W 2 dx 



(31) 



Using Eq. Q we get as the exact result the SWKB ex- 
pression 

I v / 2m(E - W 2 )dx =2irhn, n = 0, 1, 2, 3, • • • (32) 

which yields the quantum spectrum of V\ (x) . 

A similar derivation may be carried through for class 
II superpotentials obeying Eq. I|21|) . The starting point, 
as before, is Eq. I|25|) . and the integral to be considered 
is now of the form 



Jk 



Ve (A + BW 2 ) k (l 



cw 



VA+BW 2 



Ve 



VE-W 2 



dW . (33) 



The second bracketed term in the numerator on the RHS 
may be expanded binomially, and the odd-powered terms 
in W vanish on integration. We then have Jk = 



£ 



^ (fc _ 2n)!(2n)! 



Ve 



(A + BW 2 ) k - ,l (CW) 



Ve 



s/E-W 2 



-dW , 

(34) 

where n max = k/2 for k even, and (k — l)/2 for k odd. 
The highest power of W in the numerator is again W 2k 
and again only terms with this highest power (with co- 
efficient B k - n (j 2n ) will survive when Jk is differentiated 
k— times. Accordingly, Eq. (|25|l reduces to 



Si(E) = 2V2^l / VE - W 2 dx 

J Xi s 

( « (2fc-l)!! ^ 



k\ B k ~ n C 2n 



J (fc + l)(2fe)H ^ (fe-2n)!(2n)! 



35) 
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FIG. 1: Numerical evaluation of the trace formula (1 1 1't for 
the infinite square well where Fi(E) = (1 + E / E ) 1/2 - 1. 
In the figure, E is plotted in units of Eq. To ensure uniform 
lineshapes, correct degeneracies, and strict numerical conver- 
gence, we have used the usual prescription used in numerical 
semiclassics (see, for example, Section 5.5 of Ref. |8() which 
is to convolve the trace formula with a Gaussian of width a. 
For this particular calculation, we have truncated the sum at 
kmax = 10 4 while prescribing a = 0.05. 

The main results given earlier by Eqs. (JSTJ E3> remain 
valid. 

The summation in eq. (|35|) can be done similarly to 
that in Eq. (|3U[) . The inner summation provides the 
mean of (B ± Cy/B) k . Then we find 

V = 77— t 1 ~ V 1 - z +\ + 7T - 1 1 - V 1 - z -\ where 

ZZ-\- ZZ— 

z± = B±CVB. (36) 

These results (|3(JI I36[) are a simple demonstration of the 
relation between WKB and SWKB, which Barclay 4] 
approached in a different manner. 

It may now be instructive to illustrate our results with 
a few examples: 

1) Infinite Square-well. In this example, W(x) = 
— Hit / (^/ 2mL)cot(iTx / L) . It belongs to class I with A = 
h\ 2 /(2mL 2 ) = E , B = 1, and C = 0. The quan- 
tum spectrum of V± is given by f(n) = n(n + 2)Eq, with 
n = 0,1,2,.... Then Fi(E) = (1 + E/E ) 1/2 - 1. A 
careful numerical evaluation of the trace formula i|12|) 
with this Fi(E) reproduces the quantum spectrum (see 
Fig. [IJ. It is also easy to check Eq. I|31[) by evaluating 
the action integral of W 2 (x) analytically, and Eq. 11' I'D 
using Eq. (gOJ (r) = 1). 

2) 3-dimensional harmonic oscillator in the / partial 
wave. In this example W(r) = ^/2mujr/2 — %/ (y/2m)(l + 
l)/r. It belongs to class II with A = huj, B = 1/(21 + 2) 



and C = —\f~B. The quantum spectrum, measured from 
the lowest state in a fixed partial wave is f(n) = 2nhuj, 
so F(E) = E/(2hw). Again, Eq. ma y be checked 

explicitly. 

To verify Eq. (J2U, we find from Eq. JSBJl that 



V = 5 + 7^+2-V / ^+l)] ^ 

in this example. The first 1/2 represents the usual 
half-integer quantisation in LOWKB, while the terms in 
square brackets arise from the sum of order h 2 and higher 
corrections. As discussed in detail by Seetharaman 
and Barclay Q they can be removed by adopting the 
Langer prescription. We have also checked other exam- 
ples analytically. 

In conclusion, we have given a new proof that lowest 
order SWKB quantisation is exact, starting from peri- 
odic orbit theory, rather than by examining the higher 
order WKB terms. The key ingredients have been an 
invertible algebraic expression for the energy spectrum, 
and Barclay and Maxwell's 0, 0] insight about shape 
invariant potentials. 

R.K.B and D.W.L.S are grateful to NSERC for con- 
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